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Abstract 


e Numerical analysis of truncated periodic structures by utilizing conventional method 
of moments (MoM) technique is challenging because the number of unknowns 
involved is often very large; hence, the computational cost is high. This letter 
presents a novel approach, based on the dipole moment (DM) method, for analyzing 
electromagnetic scattering from truncated periodic structures whose geometrical 
and/or material properties vary locally from those of a strictly doubly infinite periodic 
structure. The main advantage of using the DM method is that it provides convenient 
closed-form expressions for the scattered fields, facilitating the matrix computation in 
the context of the MoM much more efficiently than when the conventional Rao- 
Wilton—Glisson (RWG) type of basis functions are used for the same computation. 
Numerical results obtained by using the proposed approach are compared with those 
derived from a commercial MoM software package, and good agreement is found. 
The time advantage is clearly evident from this comparison. 
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Introduction 


: Periodic structures are widely used in several applications such 
as electromagnetic bandgap (EBG) surfaces, frequency 
selective surfaces (FSSs), phased arrays, to name a few 
applications. 


e Solution of problems involving periodic structures is 
computationally expensive when conventional Method of 
Moments (MoM) technique is utilized since the number of 
unknowns is often prohibitively large. 


: To circumvent this problem, it is common to work with doubly- 
infinite periodic structures and utilize the periodic boundary 
condition (PBC) for the unit cell. 

- This assumption ignores the truncation effects and introduces 
inaccuracies in the numerical solution to the truncated periodic 
problem. 


Research Motivation 


e Solution of problems involving periodic structures is 
computationally expensive when conventional Method of 
Moments (MoM) solver ts utilized, since the number of 
unknowns Is often prohibitively large. 

e- [he proposed dipole moment approach reduces the 
computational cost by a large factor maintaining the 
accuracy of the solution results. 





Novel dipole-moment-based approach to handle 


quasi-periodic structures 


Induced current distribution in the unit-cell element is 
computed using a commercial MoM solver (FEKO). 
Unit-cell element is operating in free-space or infinite 
periodic environment. 


Far-field pattern P for the unit-cell element is 
computed using the above current distribution. 


The unit-cell element is represented by x-. y-. and z-directed 
components of the dipole moments with piecewise-sine 
current distributions, and patterns P,, P,, and P, respectively. 


The weight coefficients (Ww). w, and w3) of the three DM 

components are computed next. by matching the far fields 

produced by the unit-cell element to those produced by the 

dipole moment representations of the same using 
<P?,P> _ <PEP> 0 €PLP- 


Wi sper’ “2 ~ Serres’ 3 ^ RES 





The solution vectors of the DMs are now derived using 
the computed MoM matrix and the right hand side 
(RHS) vector. and are used to compute the fields 
scattered by the truncated structure. 


Finally. the MoM matrix elements for the problem at 
hand. namely the truncated periodic structure, are 
computed by evaluating the dot product of the closed- 
form scattered fields from the source dipole moment(s) 
and the test dipole moment(s) as illustrated in Fig. 2. 
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Matrix calculation using the proposed dipole moment approach. 





Numerical results 





Truncated periodic structure 
7x7 tilted circular loop array 


: Problem definition 
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: Center-to-center distance between the adjacent loops along the y- 
and z- directions = A/4 = 18.75 mm. 


Truncated periodic structure 
7x7 tilted circular loop array — Near fields comparison 
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E. field comparison along observation line 1 
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Truncated periodic structure 
7x7 tilted circular loop array — Near fields comparison 
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E, field comparison along observation line 1 
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Truncated periodic structure 
7x7 tilted circular loop array — Near fields comparison 


magnitude phase 











— FEKO — FEKO ; $- Interface» 

Interface size: 22.5 x 22.5 cm o 375m ; 
121 x 121 interface nodes. Qo 375cm l 
Plane wave 7 n 


Ei = 1 V/m o 
Theta — 45 degree, phi = 45 degree X — O IN 
md e 
f=4 GHz O i 
eo Loops are symmetrically 


Ez phase (degree) 


placed about the origin 
(0,0,0), and are tilted by 45 
degrees about the x-axis. 


Interface scattered fields obtained using 
DM approach and commercial MoM solver. 
az z;* Loop diameter: 5.625 mm 
( —.— Observation line | 
corresponds to z—0. 


Interface is parallel to yz plane at a 
distance of 7.5 mm from the center of 
each loop along the x-direction. 





0.05 0.1 


-0.1 -0.05 0.05 0.1 -0.1 -0.05 


0 
yo(m) 


0 
yo(m) 


E, field comparison along observation line 1 


The efficiency factor = ratio of the 
computational times of the MoM matrices 
generated by using the DM approach and a 
commercial MoM solver, is found to be on 
the order of 10°. 





Truncated Periodic Structure - Cross 
Problem definition 


: 81 x 51 truncated periodic structure with unit element as 


en 


ClOSS. 


2 
. f = 4 GHz, A = 75 mm. 

: Length of each arm (L) = A/10. 
: Incident field: Plane wave with normal incidence. 
: E; 21 V/m, theta 2 90 degree, phi 2 0 degree. 


: opacing between adjacent cross elements along y and z: 
dy = A/4, dz = NA4. 
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Coefficient vector magnitude 
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Cross element for 51 x 51 
array 





Distribution along 1st row 2-dimensional spatial distribution 
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( — 0? cut plane 
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Cross element for 51 x 51 
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Far field comparison - DM approach vs FEKO 


p = 45? cut plane 
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Matrix and coefficient vector 58.974 seconds 8.3281e+05 seconds 
computation time 


Efficiency factor = 1.4122e+04 





Truncated Periodic Structure — Elliptical patch 
Problem definition 
: 51 x 51 truncated periodic structure with unit element as 


elliptical patch. 


f24 GHz, A2 75 mm. W 
e Length of major axis (L) = 7.5 mm (A/10). aj 
: Width W(minor axis) = 4 mm. 


- Elliptical patch has a slot of length 4 mm, and is rotated about z 
axis by 45°. 


e Incident field: Plane wave. 
- E,=1 V/m, theta = 90 degree, phi = 0 degree. 


e Spacing between adjacent patch elements along y and z: dy = 
A4, dz = N/4. 
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Far field comparison - DM approach vs FEKO 
p = 45? cut plane 
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Elliptical patch 
element for 51 x 51 
array 


lime resource comparison 
DM approach vs FEKO 





Elliptical patch 
element for 51 x 51 
array 


a Ee a e 


Matrix and coefficient 58.98 seconds 8.5937e+05 seconds 
vector computation 
time 


Efficiency factor = 1.457e+04 


Observations and Conclusions | 
- The proposed dipole-moment-based approach provides 


accurate results for problems involving truncated and 
quasi-periodic structures. 


: The efficiency factor is found to be on the order of 10? for 
the /x/ truncated periodic structure problem with tilted 
circular loop as unit-cell element, and the factor is found 
to be on the order of 10^ for the 51x51 truncated periodic 
structures, for crosses and tilted elliptical patches as the 
unit-cell elements. 


- The efficiency factor becomes more favorable as the size 
of the truncated periodic structure is further increased. 
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